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ABSTRACT 

The high pericenter velocities (up to a few percent of light) of the S 
stars around the Galactic-center black hole suggest that general relativis- 
tic effects may be detectable through the time variation of the redshift 
during pericenter passage. Previous work has computed post-Newtonian 
perturbations to the stellar orbits. We study the additional redshift ef- 
fects due to perturbations of the light path (what one may call "post- 
Minkowskian" effects), a calculation that can be elegantly formulated as 
a boundary-value problem. The post-Newtonian and post-Minkowskian 
redshift effects are comparable: both are 0(/3 3 ) and amount to a few 
km s _1 at pericenter for the star S2. On the other hand, the post- 
Minkowskian redshift contribution of spin is 0{f3 5 ) and much smaller 
than the 0(/9 4 ) post-Newtonian effect, which would be ~ 0.1km s _1 for 
S2. 



1. Introduction 

Over the past decade, dozens of fast-moving stars orbiting a large compact 
mass (thou ght to be a supermassive black hole wit h mass ~ 4.4 • 1O 6 M ) have been 



discovered (IGhez et al.ll2008l ; iGillessen et al.l 120091 ) . The highly eccentric orbits and 
the low pericenter distances of these stars (~ 3 x 10 3 of the gravitational radius in the 
case of S2) provide a lucrative testing ground for general relativistic perturbations 
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to Keplerian orbits. High resolution spectral and astrometric measurements of such 
stars would also aid in the modeling of the mass distribution in the Galactic center. 
Other consequences of general relativity, the possible form of the metric itself (and 
the corresponding theories to which the metric is a solution) could be studied. 

The prospect for measuring general relativistic pericenter precession, an 0(/3 2 ) 
effect (where (3 is the pericenter velocity in light units) fo r the Galactic-center stars 
is widely appreciated (see, for example, IJaroszynskil Il998t Fragile fc Mattewsl |2000| ) 
and is anticip ated to be observable b y interferometric instruments currently under 
development (jEisenhauer et al.l 120091 ) . Should stars further in be detected in the 
future, precession e ffects at 0(/3 4 ) would become measurable, enabling tests of no- 
hair theorems (jWilll 120081 ) . 



Relativity also perturbs the kinematics of a star. IZucker et al.l ( 120061 ) drew at- 
tention to the rapid velocity changes that stars undergo around pericenter passage, 
and argued that the 0(/3 2 ) effect of tim e dilation in the s t ar's fr ame would be measur- 
able as a perturbation of the redshift. iKannan fc Sahal (120081 ) calculated the 0(/9 3 ) 
kinematic contribution of space curvature and 0(/3 4 ) effect of black hole spin g j, 
suggesting th at even the latter m ay become measurable with future interferometric 



instruments. IPreto fc Sahal (120091 ) presented a new orbit integration method in the 
presence of space curvature, spin, and Newtonian Galactic perturbations. 

In this paper, we extend previous work to include the redshift contributions that 
come from the effect of space curvature and spin on the light path between the star 
and the observer. Specifically, we will compute the redshift of a moving point source 
in a weak-field approximation of a Kerr metric, identify the various contributions, 
and investigate the scaling of these signals with orbital size. 

We formulate the problem as two photons emitted from a stellar orbit, an in- 
finitesimal proper time interval At apart. Both photons hit the observer, but with 
a difference At in arrival time. The redshift z is then 



At 



'1 + z)At. 



(1) 



The physic al process is the same as in the calculation of the spectrum of an accretion 
disk (e.g., iMiiller fc Camenzindl 120041 ) . only the computational strategy needed is 
different. In the accretion- disk case, photons are shot backwards from the observer 
in a range of directions, but in the general direction of the extended source (the disk). 
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By considering the collision events of each ray with the source surface, the quantities 
of interest may be read off. In the stellar-source case, however, it is necessary to 
solve for the initial direction of a photon so that it will reach the observer. Hence 
we have a boundary-value problem. 

Figure [1] is a schematic illustration of our method. Photons are emitted from 
nearby spacetime points on the stellar orbits and travel to the observer. In the left- 
hand picture, the star emits "Minkowski photons" which feel no space curvature and 
travel in straight lines; redshift depends only on the velocity and time dilation of 
the star. This is in effect the approximation used in previous work. In the middle 
picture, the star emits "Schwarzschild photons" which feel space curvature. In the 
right-hand picture, the star emits "frame dragging photons" which feel spin as well 
as space curvature. 

Below, § 12. ll details the problem to be solved and the method used for calculation 
of the redshift. The Matlab scripts implementing our algorithm are available as an 
online supplement. Then § 12.21 presents the black-hole model and associated metrics 
which we use in our approach, and § 12.41 derives how the various effects scale with 
orbit size. We apply our algorithm to the star S2, and detail the results in Section [3j 



2. The redshift-calculation method 

2.1. A boundary- value problem 

In order to calculate the redshift of a moving star as observed by a fixed observer, 
we need to solve the geodesic equations for both the star and for photons. Geodesic 
equations are commonly expressed in terms of the Lagrangian, 

with dots denoting derivatives with respect to the affine parameter. But an equivalent 
formulation exists in terms of a Hamiltonian 



(3) 



-4- 



We will follow the latter in this work. Numerically 7i = C. We will in fact employ 
two Hamiltoniana^ "Hstar and "H nu11 , which are really two different approximations for 
the same Hamiltonian. Numerically "H nu11 = of course. We will write A for the 
affine parameter of the star, and a for that of a photon. 

Consider two photon trajectories, emitted at two points on the star's orbit AA 
apart in the affine parameter, with both photons terminating at the observer and 
arriving at time At apart in the observer's frame. Since dr = ^y\g fMU dx^dx u \, we are 
able to express the proper time between the emission events in terms of the affine 
parameter as 

At = AAv/|2H st ar| (4) 

where ?4tar is evaluated at either point of emission. Comparing with the definition ([Q) 
of the redshift, we have 

At 

z = . - 1. 5 

A\^/\m^\ 

We now need to compute At for a given AA. To do this, we begin by calcu- 
lating the orbit of the star for some chosen initial conditions, by solving Hamilton's 
equations for "H sta r with A as the independent variable. The temporal component of 
the generalized momentum is set at pt = — 1. This amounts to choosing the units 
for A such that dt = dX outside of gravitational fields. We then choose a point on 
the star's orbit whose observed redshift we wish to calculate, and from this point we 
seek a photon that will reach the observer. 

Consider a function $ J , which effectively shoots a photon by integrating Hamil- 
ton's equations for "H nu11 with given initial conditions at affine parameter a = and 
returns the 3-position at a = 1. We write 

&(t,x i ,p i )=x i \ 0=1 (6) 

where i, j = 1,2,3 and the initial p t is chosen such that H nu11 = 0. We pass the 



*In an effort to keep the distinctions between effects on the star orbit and those on the light path 
as perspicuous as possible, we adopt subscripts for orbital effects and superscripts for light-path 
effects. This notation may be found on Hamiltonians T~L and redshifts z, and has nothing to do with 
covariant / contravariant indices. 
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function & to a root-finder, and solve for the root of 

/ (ft) = & - xi hs (7) 

by varying the initial 3- momentum ft. Naturally, we may not adjust the x\ as we 
are interested in a specific point on the star's orbit. 

The root-finding algorithm requires a set of initial guesses for the initial ft of 
the photon. On this account, we shoot an initial-guess photon from the star position 
in the direction of the observer, that is, we evaluate & with trivial initial conditions. 
These initial-guess values for the ft shoot in the direction of the observer ignoring 
curvature. However, because the spacetime is indeed curved, this photon will not hit 
the observer. It serves only to start the root-finding algorithm. 

Once the root-finder reports a solution within specified tolerance level, we move 
the star a very short distance (AA in affine parameter) along its orbit, and repeat 
the above procedure, now solving for the sought-after trajectory at the star's new 
position. The difference in arrival times of the two photons is At, which we insert 
into (jSJ) to evalate the redshift. 

Thus far we have solved for two photon trajectories. This computation has 
enabled us to calculate the redshift at a chosen point on the orbit. We may then 
repeat this process at further points along the orbit, garnering results as much as 
the required resolution demands. 

2.2. Post-Newtonian and post-Minkowskian approximations 

The spacetime outside a spinning black hole is described by the Kerr metric. 
Since the Galactic-center stars are far from the horizon, approximations valid only 
at large r can be used to study general relativistic effects. Accordingly, we first 
derive two perturbative Hamiltonians, a post-Newtonian "H sta r for stars, and a post- 
MinkowskiarJl] approximation W nul1 for photons. Different physical effects come into 



^The term 'post-Minkowskian' is often used as a synonym for 'weak-field metric'. We are using 
it, however, to refer specifically to light paths that deviate slightly from special relativity. 
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play at different orders of the perturbation parametei@ e, and we will show these 
numerically by toggling different terms on and off. 

Taking the Kerr metric in Boyer-Lindquist coordinates with geometric units 
GM = c = 1 (leaving us with a unit of length equal to the gravitational radius 
GM/c 2 ~ 5 x 10 6 km for the Galactic-center black hole) we have the full Hamiltonian 

_ (r 2 + g 2 ) 2 - S 2 Asin 2 fl A 1 

A - s 2 sin 2 9 o 2sr . . 

"WW + va w * (8) 

(9) 

where 

A = r 2 - 2r + s 2 and p 2 = r 2 + s 2 cos 2 9. (10) 
and s denotes the black hole spin parameter. 

Let us first consider the dynamics of the star. Sufficiently far from the black 
hole, the post-Newtonian approximation 

v 2 ~l/r (11) 

applies. If we choose v ~ (e) then by ffTTj) r is (e~ 2 ). Correspondingly, we force 
the velocity terms p r , pe/r and p^/r to be (e). Making the following replacements 
in (JED 

r -)■ e~ 2 r, p r -)• ep r , p e -)■ e~ l p B and p -)■ e~% (12) 
and keeping terms to (e 5 ), we obtain 



2 \ 2 2r 2 2r 2 sin 2 9 r J \ r 2 r J r 3 
We can abbreviate ( fl~3l) as 

%star = "^static + ^ "^Kop + ^ "%Schw + £ 5 "%FD- (14) 



^In this paper e, e 2 and so on are just labels to keep track of different orders. Numerically e = 1. 
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Here ?4tatic produces motionless geodesies, "Hxep gives the Keplerian phenomenology, 
"Hschw is the weak-field Schwarzschild contribution that produces pericenter preces- 
sion, while the angular-temporal term "Hfd produces the Lens-Thirring effect or 
frame dragging. 

Continuing now to photon trajectories, we remark that by the equivalence prin- 
ciple, the full Hamiltonian is exactly the same for photons and stars. However, the 
same terms can have different orders in the two regimes, prompting approximation 
"H nu11 . In particular, the approximation (fTTj) obviously does not hold for photons, 
which have v 2 = 1, implying the scalings 

r -> e~ 2 r, p e -> e~ 2 p e and p^ -> e~ 2 pj, (15) 

with p r being 0(1). Expanding as before, we obtain 

a/null = _P± vl Ve_ V\ (Pt Pr\ 2 | 

2 2 2r 2 2r 2 sin 2 fl V r ^ / 

2p t 2 2sm^ s 2 sin 2 ^ 2 s 2 cos 2 # 2 \ 4 

+ — 5 Pr + „ 4 Pe e ( 16 ) 



r 



2 1 r 3 2r 2 2r 4 



- a different selection of terms compared to the post-Newtonian case. There is no 
term at 0(e 5 ). We may abbreviate (|T6|) as 

U nul1 = U Mink + e 2 -H Slo + e 4 (^ Snlo + U FD + -H torq ) . (17) 

At zeroth order we have special relativistic or Minkowski photons. The leading 
order Schwarzschild effect "H Slo gives the gravitational deflection of light. At (e 4 ), 
however, there are three distinct effects: first "H Snlo gives a next-to-leading order 
correction to the Schwarzschild effect, the off-diagonal term "H FD gives frame dragging 
again but for photon trajectories, while 7/ torq provides a torque proportional to s 2 . 



2.3. Pseudo-cartesian coordinates 

The spatial Boyer-Lindquist coordinates r, 9, are convenient for computing 
stellar orbits, but not well suited for photon paths. The photon paths are nearly 
straight lines, but since the observer is much further from the black hole than the 
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source, tiny variations in 9 and <j) at the observer imply large distances. As a result, 
both the integrator and the root-finder become susceptible to roundoff error. 

To cure the problem, we change to pseudo cartesian coordinates. They are not 
purely cartesian, as the surface x 2 + y 2 + z 2 is not spherical. 



x = r sin 9 cos 0, y = r sin 9 cos <f>, z = r cos 9 



(18) 



The corresponding momenta are readily derived by completing the canonical trans- 
formation, leading to the usual relations 



(19) 



The form of H nul1 changes accordingly. The Minkowski part bears the familiar form 

,2 



2 

1/Mink _ Pt , P 



The leading-order Schwarzschild terms are 

■^Slo _ Pt 

and the associated 0(e 4 ) term retains its previous form of 



(x ■ P) 1 



2p 2 



Also at 0(e 4 ) we have the frame-dragging term 

2sp t 



and the torquing terms 



^torq 



S 

2^2 



1 - 



n FD = 



xp 



(x x p) 2 , 



s 

2? 



P 2 - 



x P 



y _ lx x p) 

J r 2 _ 



z* 



(20) 



(21) 



(22) 



(23) 



(24) 



The torquing terms are the most difficult to integrate numerically. This is because 
while the terms themselves are at (e 4 ), they involve quotients of particularly high 
powers. Taking derivatives of "H torq further bloats these bottom and top heavy frac- 
tions, and provokes roundoff errors. We argue below that the torquing terms are 
in any case unimportant for the known Galactic-center stars. Hence we omit these 
terms in our numerical work. 
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2.4. Scaling properties of redshift contributions 

We can infer the scaling with orbital size a of the redshift contributions Az of 
the various perturbative terms in "H s tar with the following deliberation. Consider a 
perturbative term 

A^ star ~ a- n (25) 

Since "H sta r is constant along the orbit, any variation in A"H star has to be balanced by 
a variation in the unperturbed Hamiltonian. Since the latter scales as "H s tar ~ l/o, 
we have AT-t sta r ~ Aa/a 2 , giving 

Aa ~ a 2 - n (26) 

Furthermore, since the orbital velocity scales as v ~ a" 1 / 2 , we have At> ~ ar 3//2 Aa 
and hence the redshift signal 

Az (A# star ) ~ ah- n . (27) 



A similar argument can be made for the perturbative terms A"H nu11 . In this case 
we compare photons emitted from the same point, only with different Hamiltonians. 
While the redshift signal from the previous case necessitated our consideration of the 
stellar velocity only, in analyzing the gravitational redshift signal, we are naturally 
interested in the time difference At due to AH. Let the perturbation be 

A7* nuU ~r- n . (28) 

Since the light travel time is an integrated quantity, we expect the change in the 
light travel time to scale as 

At trav ~ r 1 -" (29) 

The change At trav in light travel time must not be confused with the difference At 
in arrival time of two photons. For the latter quantity we get At ~ Ar/r n and since 
Ar between two photons is ~ v ~ a -1 / 2 we derive 



Az (AU nul1 ) ~a~^ n 



(30) 
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Using ([27)1 and ([30]) and recalling that each power of e in the Hamiltonians 
represents a scaling factor of r -1 / 2 we can read off the following scalings. 



%Schw 
%FD 

■^torq 



FD 





/ 


a" 


-3/2 \ 






a' 


-2 


oc 




a' 


-3/2 




V 


a 


-5/2 



(31) 



3. Application to S2-like orbits 

We select S2 for study, since S2 has the shortest orbit and one of the 

highest peri center velocities of all the known stars orbiting Sagittarius A*, and hence 
provides us with perhaps the best opportunity to observe general-relativistic effects. 

Figure [2] shows a redshift calculation f or a star with S2's orbital parameters. 



These are taken from iGillessen et al.l ( 120091 ) . The gravitational radius is taken as 
5 x 10 6 km. For defmiteness, we take the spin to be maxim al, and pointing towards 



Galactic North. Disc seismology models ( iKato et al.ll2009l ) put the spin at s ~ 0.44. 
The direction however, remains unknown. 



All the contributions to "H star in ( TT41) are included for the orbit calculation. The 
photon trajectories include all contributions to "H nu11 in (fT7j) except for "H torq . 

Naturally we would like to compute the redshift contributions of the various 
relativistic terms, and verify that they follow the expected scalings (I3T|) . In order 
to do this, we examine the differences between redshifts computed from different 
post-Newtonian and post-Minkowskian cases. This allows us to isolate the effects of 
Hschw, H FD , % Slo and H Snlo + H FD , as follows. 

1. To isolate "Hschw we compute the redshift difference — ^Kep k - By we 
mean that the star is followed using terms in "H sta r up to "Hschw and the photons 
are followed using 7-L nuU up to 'H Mmk . The same naming convention applies to 
z Kep k ari d to other expressions of this type below. 

Figure [3] shows the redshift difference, calculated for three orbits going from 
apocenter to apocenter. One orbit has the parameters of S2; the two others 
have a = 2a§2 and |as2 with the other orbital parameters being the same. The 
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redshift difference increases till pericenter and then declines somewhat, but 
not to its previous apocentric value, because the relativistic orbit experiences 
prograde Schwarzschild precession whereas the Keplerian orbit does not, and 
the resulting phase change in the orbit gives an increasing contribution to the 
redshift. 

For S2 parameters, the maximum redshift contribution of T^schw is found to to 
be ~ 7 km s _1 . For the other two stars, upon rescaling the redshift differences 
by (a/as 2 ) -3 / 2 and the orbital time also by (a/cts2)~ 3//2 , the results can be 
overlaid almost perfectly on those of the S2-like star. 

2. To isolate Hfd we then compute the redshift difference z^^ k — z$^£. For an 
S2-like orbit the signal is around 0.1km s _1 at pericenter, and as Figure S] 
shows, the signal scales as a -2 . 

3. To isolate H Slo we compute Zgchw — z<£^ and illustrate this difference in Fig- 
ure |5j There is no precession-related redshift effect involved, because the pho- 
ton types being compared refer to the same stellar orbits. The signal scales as 
a~ 3 / 2 and for the S2-like orbit the maximum is ~ 2 km s _1 . 

We see that the Schwarzschild terms in the stellar orbit and in the light path 
give comparable contributions to the redshift. To detect the Schwarzschild 
effect, it is necessary to take both into account. 

4. Finally, we isolate "H Snlo +"H fd by computing — z^ . Figure E] verifies the 
expected a _5//2 scaling and shows that the maximum signal for S2 parameters 
is ~ 10~ 2 km s _1 . Thus we see that the frame- dragging on photons is much 
smaller than on stars. Similarly, % Snlo makes a much smaller contribution than 
Hfd- We expect the contribution of % torq would be similarly small, though we 
have not calculated it. 

We see that in order to measure the leading-order frame- dragging effect on 
Galactic-center stars, it is sufficient to consider Schwarzschild photons. 



Of course, the computation method for weak redshift signals contains numerical 
errors, especially for higher-order effects being evaluated at large distances from the 
black hole. The numerical noise in our implementation drowns out the the leading- 
order Schwarzschild effect for scaled S2-like orbits with a ~ 3 • 10 8 — at which point 
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the pericenter redshift signal is ~ 5 • 10~ 6 km s _1 . Numerical noise would overwhelm 
the frame-dragging signal for scaled S2-like orbits with at a « 1.2 • 10 5 . This being 
an order of magnitude larger than 052, and so we are in good shape to calculate 
the frame- dragging redshift contribution. For redshift signal contributions beyond 
those of frame-dragging, the numerical noise in our Matlab implementation of the 
algorithm is intolerable for 052. Were we calculating these effects for stars closer 
to the SBH, the higher-order signals would be stronger, and therefore less prone 
to round-off. In summary, for the known Galactic-center stars, the numerical noise 
in calculating relativistic redshifts will be well within the observational errors, even 
with the next generation spectral instrumentation. 



4. Summary and Outlook 

Some stars in orbit around the Galactic-center black hole reach velocities of a 
few percent light at pericenter, and the time-varying redshift of these stars during 
pericenter passage has small but distinctive perturbations from general relativity. 

The redshift is dominated by the line of sight velocity, which for the star S2 
reaches v ~ 5 x 10 3 km s" 1 at pericenter. The leading perturbations are from time 
dilation because (a) the star is moving, and (b) because it is in a potential well. Both 
of these make 0(/3 2 ) contributions to t he redshift (where y is the stellar velocity in 
light units) and are well understood (IZucker et al.l 120061 ). In this paper we have 



calculated additional perturbations from general relativity, which are the following. 

1. The weak- field Schwarzschild effect on the stellar orbit, which contributes to 
redshift at 0(/3 3 ). For S2 it is ~ 7km s" 1 ! 

2. The frame-dragging effect of black hole spin on the stellar orbit, which perturbs 
the redshift at 0(f3 4 ). For S2 it would be ~ 10 _1 km s _1 for maximal spin. 



§Note that these effects are not measurable separately, only the total redshift is. Hence the 
values like 7 km s^" 1 depend on the choice of reference orbit and phase, and can change accordingly. 
Nevertheless, the stated numbers give an idea of the observational precision required. 
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3. The weak- field Schwarzschild effect on the light travelling from the star to us, 
which gives a redshift perturbation at 0(/3 3 ). For S2 it is ~ 2km s _1 . 

4. Frame dragging plus next-order Schwarzschild perturbation of the photon 
paths. These contribute at 0(/3 5 ) to the redshift, and we estimate these as 
~ lO^kms- 1 for S2. 



Of these, the first two are orbital effects an d have been considered in previous 
work ( iKannan fc Sahall2008l ; iPreto fc Sahall2009l ). The last two are light-path effects, 
known about but not previously computed in the context of Galactic-center stars. 

Clearly, in order to measure the Schwarzschild effects via the redshift, the effects 
of general relativity on both the stellar orbit and the light path must be computed, 
as they are of the same order. On the other hand, to leading order, frame dragging 
needs to be considered only in the orbit and can be neglected in the light path. 

In order to test for the the presence and form of the NLO and NNLO terms in 
the metric, the calculated redshift curve must be fitted to the spectral data via a 
range of parameters. These include the orbital parameters and the black hole mass. 
In keeping these parameters variable, we can expect a requirement for spectroscopic 
accuracy less than the signal sizes themselves. Bear in mind however, that this 
depends on the number of data points. The inclusion of astrometric data to the 
fitting procedure will help relax the accuracy bound. 

Observationally, the Galactic-center stars are of course observable only in in- 
frared. For the most massive stars, including S2, the Brackett-7 line at 2.16 /im 
is the most prominent available spectral fea ture, and has an intrinsic dispersion of 



100km s 1 (see e.g., iMartins et al.l 120081 ) . For low-mass stars, the edges of the 



CO molecular bands (the so-called CO band heads) are excellent sharp features for 
redshift measurements. 

The SINFONI spectrograph, an instrument at the VLT, has a spectral resolution 
of 75 km s -1 . Redshift e rrors are currently estimated at 10s of km s _1 (see Section 4.1 
of iGillessen et al.l 120091 ) but expected to improve. Measurements by this instrument 
during S2's next pericenter passage in 2016 could suffice to provide data from which 
the Schwarzschild signals could be extracted. 



Spectral measurements seeking the spin-dependent signals are far beyond the ca- 
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pabilities of existing infrared spectrographs. On the other hand, re cent developments 



such as laser-comb spectrographs (see e.g.. ISteinmetz et al.ll2008l ) suggest optimism 



that the spectral resolutions of next-generation instruments will prove adequate. 
Meanwhile, there are some theoretical issues that require further research. 

• Given a model of gravity which is metric, and an associated energy-momentum 
distribution, for which we have a metric, we are able to calculate the general 
relativistic redshift as observed by the Earth. Should we wish to probe the 
agreement of measured redshift contributions with a model, an effective way 
of working backwards needs to be formulated. Given redshift curve data of 
sufficient resolution, methods for determining the model from such must be 
developed. 

• The Kerr Black Hole metric is a vacuum solution to the Einstein field equa- 
tions. The assumption of a 'clean' metric is an oversimplification. An extended 
Newtonian mass distribution in the galactic center is anticipated. Accreting 
material, gas, and a possible accumulation of dark matter in the center could 
play a significant role in the dynamics. Such mass distributions need to be 
included in the metric. 

• The possibility of a non-Einsteinian black hole metric should not be disre- 
garded. Alternative theories of gravity possess black h ole solution s whose phe- 



nomenology differs potentially already at low order (jWilll 119931 ) from weak- 
field Einstein gravity. Suitable measurements, combined with the concession 
of extended mass distributions would allow for the identification of potentially 
revealing redshift contributions. Conclusive results of such studies would likely 
require spectral measurements with a resolution beyond that offered by present- 
day generation instruments. 



We thank S. Gillessen and the referee for comments. 
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Fig. 1. — Schematic illustration of the method. Pairs of photons are emitted at 
slightly different proper times along the orbit, in precisely the right direction to 
reach the observer. Finding these photons is a boundary-value problem, and once 
found, each photon pair allows us to calculate the redshift at that point on the orbit 
by evaluating ([5]). At the left of the figure we have Minkowski photons, which move 
in straight lines. In the middle we have Schwarzschild photons, which are lensed. On 
the right, we have frame-dragged photons. The time difference between the emission 
of each photon in a pair has been exaggerated here for visual clarity. Note that it is 
only the star's unrealistic proximity to the black hole that allows for such a visible 
depiction of the different effects. 
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Fig. 2. — A redshift calculation for an S2-like star over two orbits. The left panel 
shows the time each photon takes to travel from the star to the observer. The star 
begins at apocenter, which happens to be closer to the observer than the pericenter. 
Therefore the light travel time naturally increases as the star moves towards its 
pericenter. The right panel shows the redshift. The peaks occur during pericenter 
passage due to the high pericenter passage velocities. The calculation includes post- 
Newtonian and post-Minkowskian effects, but these are too small to see in this figure. 



- 18 - 



CO 

E 

CD 
O 
C 
CD 

CD 



CO 
T3 

CD 
DC 




5000 10000 1000 2000 1000 2000 
Emitted Time (days) Shifted Emitted Time (days)Shifted Emitted Time (days) 



Fig. 3. — Redshift difference - z^ k showing the contribution of %schw The 
left panel shows the redshift difference for three orbits, one with the parameters of 
S2, and the other two with a doubled or halved. In the middle panel, the horizontal 
scale is stretched (a/a S 2)~ 3//2 - In the right panel the redshift difference is scaled by 

(a/a S2 )" 3 / 2 . 
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Fig. 4. — Redshift difference 2^ nk — ^schw shewing the contribution of "Hfd- The 
scheme follows Figure [3] except that the redshift difference is scaled by (a/a^)^ 2 - 
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Fig. 5. — Redshift difference — ^sSw showing the contribution of "H Slo . The 
scheme follows Figure [3j the redshift difference being scaled again (a/as2)~ 3 ^ 2 - 
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Fig. 6. — Redshift difference Zp® — z F ^ showing the contribution of "H Snlo + "H FD . 
The scheme follows Figures [3H5J with the redshift difference scaled by (a/as2)~ 5 ^ 2 - 
Here we have taken the spin as maximal, s = 1. This signal is exactly proportional 
to s. 



